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In this paper, we investigate the dynamics of synchronous totally asymmetric exclusion processes 
(TASEPs) on lattices with a multi-input-single-output (MISO) junction, which consists of m sub- 
] chains for the input and one main chain for the output. An MISO junction is a type of complex 

geometry, which is relevant to many biological processes as well as vehicular and pedestrian traffic 
■ flow. A mean-field approach is developed to deal with the junction that connects the subchains 

| and the main chain. Theoretical calculations for stationary particle currents, density profiles and a 

phase diagram have been obtained. It is found that the phase boundary moves toward the left in 
the phase diagram with the increase of the number of subchains. The non-equilibrium stationary 
states, stationary-state phases and phase boundary are determined by the boundary conditions of 
the system as well as by the number of subchains. The density profiles obtained from computer 
simulations show a very good agreement with our theoretical analysis. 
PACS numbers: 05.70.Ln, 02.50.Ey, 05.60.Cd 
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Non-equilibrium transport phenomena have attracted much attention of physicists since physical principles under- 
■ lying these phenomena could be revealed in terms of phases and phase transitions [3, H| • Totally asymmetric simple 
exclusion processes (TASEPs) serve as a basic model for non-equilibrium systems and have been widely applied in the 
r/2 . study of transport phenomena in chemistry, physics and biology, for example, particle diffusion through membrane 

channels the kinetics of synthesis of proteins 0], polymer dynamics in dense media Q, gel electrophoresis 
43 intracellular transport of motor proteins moving along cytoskeletal filaments , vehicular traffic [3, H and ant traffic 

TASEPs are non-equilibrium one-dimensional lattice models in which particles move along one direction with hard- 
core interactions. The exact solution of TASEPs has been obtained by using the matrix product ansatz (MPA) [10| and 



the Bethe ansatz 11 L respectively. Recently, some extensions of TASEPs focus on coupling with Langmuir kinetics 
0, 0, ES H HI til El HI [111 and/or latticegeometries [H, SI SI HE SI S3, HI , as well as bottleneck-induced 
transport phenomena such as in [29l 1301 ] . Ref. [291 ] investigated diffusive compartments as bottlenecks in a driven 
transport. They found that when a diffusive bottleneck is at the boundary of a system, the system cannot reach a 
maximal current phase; when a diffusive bottleneck is in the interior of the system, the system is dictated by the 
diffusive bottleneck and has a maximal current defined by the bottleneck. More recently, Ref. [3(| introduced a 
bottleneck phase to describe the phenomenon that the current is independent of boundary conditions. 

In this paper, we focus on one special lattice geometry - junctions, which are widely observed in many real physical 
systems. Those junctions are formed for various reasons. For example, (i) variation of the number of protofilaments 
on a microtubule in vitro [3lj : (ii) transport of vesicles in a branching axon or dendrite [32| : (iii) merging of two or 
more roads; (iv) data through hubs (e.g., switches, routers) on the Internet. This kind of lattice geometry can cause 
congestion, e.g., in the traffic of molecular motors, vehicles or data packets. Traffic congestion of molecular motors 
could lead to some human diseases such as Alzheimers [33[ and some neurodegenerative diseases Q • Vehicular traffic 
congestion can pollute environment, increase fuel consumption. 

Inspired by this wide range of possible applications, we investigate the dynamics of synchronous (i.e., in parallel 
update) totally asymmetric exclusion processes on lattices with a multi-input-single-output (MISO) junction (see Fig. 
1). The parallel updating procedure has been typically adopted in modeling vehicular and pedestrian traffic [3, @|. 
Multiple inputs exhibit more complex interactions between particles at junction points than two inputs. In reality, it 
can be observed that several traffic lanes merge into one lane and multiple protofilaments come together to form one 
protofilament (3l| . However, they have not been understood well from the viewpoint of theoretical analysis. 

Theoretical calculations, along with a mean-field approximation, are developed. The phase diagram is presented 
and density profiles are investigated. A phenomenological domain wall theory, based on Refs. [H Hi|, is used to 
predict phase boundaries. Computer simulations are also conducted. Note that TASEPs on lattices with Y-junctions 
(e.g., two-input-single-out junctions) in a random updating procedure has been investigated in Ref. (23j . Y-junctions 
can be seen as a special case of multi-input-single-output junctions. 
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The paper is organized as follows. In Section II, we give a description of a synchronous TASEP model with an MISO 
junction, theoretical calculations as well as the mean-field approximation are developed. In Section III, we analyze 
the phase boundaries using a phenomenological domain wall theory. In Section IV, the results of our theoretical 
calculations and computer simulations arc presented. Finally, we give our conclusions in Section V. 



II. THE MODEL AND MEAN-FIELD APPROXIMATION 



An MISO junction is illustrated in Fig. 1. The system consists of m subchains for input and one main chain (chain 
m + 1) for output connected by junction points-sites N on the subchains and site N + 1 on the main chain. Each 
subchain and the main chain includes N sites. For simplification, inter-lane transitions between subchains are not 
permitted in this paper. Particles are assumed to move from the left to the right. 
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main chain m+1 



^ • • • _ • 
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(a) 



(c) 



FIG. 1: (a) Schematic diagram of a synchronous TASEP with an MISO junction. Particles move from the left to the right 
with hard-core exclusion, (b) In a subchain, the injection rate at site 1 and the ejection rate at site N are given by a and /3 e //, 
respectively, (c) In the main chain, the injection rate at site N + 1 and the ejection rate at site 2N are given by a e f / and j3, 
respectively. 

An occupation variable, Ttj, denotes the state of the ith site in the Ith subchain (£ = 1,2, ...,m) and the main 
chain {£ = m + 1). Tg.i — 1 (or re t i — 0) means that site t#,, is occupied (or empty). The system updates all particles 
in parallel by the following rules (see Fig. 1): 

• 2 = 1. (i) If t^i — 0, a particle enters the system at rate a; or (ii) If rt,\ — 1 and = 0, then the particle 
at site (£, 1) moves into site (I, 2); or (iii) If both Ti t \ = 1 and 7^2 = 1, then the particle at site (£, 1) does not 
move. 

• 2 = N. (i) If sites N of k subchains (1 < k < m) are occupied by k particles at the same time, particles have 
the same priority to hop to site N + 1. However, only one particle will enter site N + 1 at any single time step, 
providing that site N + 1 is empty; or (ii) If only one site N of the subchains is occupied by a particle, the 
particle can directly hop to site N + 1 providing that site N + 1 is empty. 

• 2 = 2N. If T m+ i,2N = 1, the particle leaves the system with rate f3. 

• l<2<A r orA r + l<2< 2A^. If = 1, the particle can move into site (£, i + 1) providing r^j+i = 0. 
Otherwise, the particle cannot move. 

Exactly solvable results of an one-dimensional synchronous TASEP have been obtained in Ref. [3?], [HI ■ We briefly 
recall these results, as the solution of our synchronous TASEP with a junction can be derived from them. There are 
three phases (low density (LD), high density (HD) and maximal current(MC)) and a transition line when a = (3. The 
MC J = 0.5 can only be reached at a = (3 = 1 [381 ] . 

• When a < j3 < 1, a low-density (LD) phase is obtained with 

a a 
1 + a p(l + a) 

where J is the system current; p is the bulk density; p\ (pn) is the particle density at the first (last) site. 
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When 1 > a > 0, a high-density (HD) phase is obtained with 



j = i-p, p = ^ ^ = 1 -^rrjy PN = p - (2) 

• When a = < 1, a transition line between LD and HD is obtained. 

• When a = (3 = 1, the maximal current (MC) is obtained and J = 0.5. 

Based on the above results, we develop exactly solvable results for TASEPs with an MISO junction. For an MISO 
junction, as the current is conserved through the system, we have: 

Ji + Ji + ■ ■ ■ + Jm — Jm+l, Ji = Ji = ■ ■ ■ = J m , mJe — J m +l < 0.5 (3) 

where (£ = 1, 2, m) is the current on the £th subchain; J m+ \ is the current of the main chain. 

Each of the m subchains of the MISO junction can be seen as a synchronous TASEP with injection rate a and 
ejection rate e ff, while the main chain can be seen as a synchronous TASEP with injection rate a e // and ejection 
rate (3. a e f f and (3 e / / can be written as 

Peff = 1 - Pn+i, a e fj = mp N . (4) 

These m subchains should have the identical phases when particles on the m subchains merge into the main chain 
with the same priority. Our computer simulations also support this prediction. Thus, the stationary state of the 
system can be obtained by combining the possible phases that exist in each of these subchains and the main chain. 
As each single chain may have three possible phases (LD, HD and MC), due to the equivalence of these subchains, 
the number of possible stationary phases of the system is equal to 3 2 = 9. In other words, a stationary state can be 
one of the following nine phases: the (LD, LD), (LD, HD), (LD, MC), (HD, LD), (HD, HD), (HD, MC), (MC, LD), 
(MC, HD), and (MC, MC) phases. 

One can see that three phases cannot exist: (MC, LD), (MC, HD) and (MC, MC). According to Eq. (3), it is 
impossible for the maximal current phase to exist in a subchain since the maximal possible current in the system is 
no more than 0.5. Therefore, the number of the possible phase combinations reduce to 6, i.e., the (LD, LD), (LD, 
HD), (LD, MC), (HD, LD), (HD, HD), (HD, MC) phases. 

• The (LD, HD) phase. The conditions for this case are as follows: 

a < e ff, a eff > (3. (5) 
From Eqs. (1) and (2), the stationary properties of this phase are given by: 

a a 

■Jl — r. — ; , -Jm+l = . , a i Pi — 



1 + a 1 + (3 1 + a 

a 1 1 ra\ 

PN = -5 77— 7, PN+1 = 1 n m , P2N = 1 ^ W 



According to Eq. (3), mJi = J m +i, we get: 



m+{m-l)(3 {) 



However, a e ff and (3 e ff are not solvable from above equations. In other words, we cannot calculate the bulk 
density through the above equations. The density will be solved through the domain wall theory in Section III. 
From Figs. 2(a) and (b), one can see that a = (3/[m + (m — 1)(3] (when < 1) corresponds to the transition 
line between the (LD, LD) phase and the (HD, HD) phase. 
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• The (LD, MC) phase. This phase corresponds to the following conditions: 

a<0eff, a eff =P = l. (8) 

According to Eqs. (1) and (3), we obtain: 

a =2^~r Jl = i' J ^- - 5 - ( fl ) 

Obviously, a = 1 /(2m - 1) < 1/2 when to > 2. That is, Eq. (8) is satisfied. Thus, the (LD, MC) phase can 
exist in the system when: 

a = 2m~~T' /? = L ^ 

Again, we cannot calculate the bulk density through above equations. The density will also be solved through 
domain wall theory in Section III. From Figs. 2(a) and (b), one can see that the (LD, MC) phase is the transition 
phase between the (LD, LD), (LD, HD), (HD, HD) and (HD, MC) phases. 

• The (LD, LD) phase. The following conditions should be satisfied: 

a < (3 eff , a e ff < (3. (11) 

According to Eq. (1), the stationary current and density are given by: 

a a e ff a 
J i — i — ; — , Jm+i — -r- , Pi — 



-I J "771+1 -i . ) fl — -i . J 

1 + a 1 + a e ff 1 + a 

a a eff a eff , 10 n 

/3 e ff{l + a) 1 + Oteff p{l + a e ff) 

Using Eqs. (3) and (4), a e ff and /3 e // are expressed as: 

ma a 1 - (m- l)a 

aeif= l-( m -lW 0eff= l + a ■ (13) 

Since a e ff < 1 and a e fj = ma/[l — (m — l)a], a < l/(2m — 1). Substituting Eq. (13) into Eq. (11), we 
obtain a < + m 2 /4 — to/2 for a < (3 e ff, and a < (3/[m + (m — l)/3] for a e // < /3. Since l/(2m — 1) < 
^/l + to 2 /4-to/2 (when to > 2) and /3/[m + (to - < l/(2m- 1) (as (3 < 1), the system is in the (LD, LD) 
phase when: 

a< /3<1. (14) 

to + (m — l)p 

The (HD, HD) phase. The conditions for this case are given by: 



a > {3 e ff, a eff > f3. (15) 



The current and density of this phase in a stationary state are: 



j _ Peff j _ 13 _ 1 _ Peff 

I + Peff 1+P a(l+/3 e //) 
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^ w = T~T« — ' P N+1 = 1 n _l m ' ^ 2Ar = TXfl- ( 16 ) 

i + Peff ot e ff{l+p) l+p 

From Eqs. (3) and (fTB]) , we obtain /3 e // = /3/[m + (m — 1)0]. Thus, the system is in the (HD, HD) phase when: 

a > / T\a ( 17 ) 

m+ (m-l)/3 

• The (HD, MC) phase. The corresponding conditions for this phase are: 

a>p e ff, a e //=/3 = l. (18) 

According to Eqs. (2-3), we obtain 

1 1 2m- 2 1 

Jl = 2W PAr+1= 2^T' ^=2^1- (19) 



Thus, the (HD, MC) phase can exist in the system when: 

a>—?—, [3 = 1. (20) 
2m — 1 

(HD, LD) phase. The conditions of existence for this phase can be written as 

a>/3 ef f,a eff <p. (21) 
The corresponding expressions for stationary current and density are 

J\ — - — —3 — , Jm+i — "j — ; 1 Pi — 1 



1 + " l l + a e //' a(l + /3 e //)' 

PN -TT]^> ^'TT^' ^ = /3(i + a e//) - (22) 



According to Eq. (3), we have 



From Eqs. (4) and (22), a e // and $»// can be rewritten as follows 



in 1 

»e/J ~ 1 I a . &ff = T~Z • ( 24 ) 
1 + fj e ff 1 + a eff 



Substituting Eq. (24) into Eq. (23), we obtain a e // = y/l + m 2 /4 + m/2 - 1 and /3 e // = \/l + m 2 /4- m/2. It can 
be seen that values of Q. e ff and f3 e ff are determined by the number of subchains m, independent of a and /3. This 
indicates that the (HD, LD) phase cannot be represented in the a — f3 plane. In other words, the (HD, LD) phase 
does not exist in the system. In fact, when the subchains are in the high density phase, particles at site ./V will hop 
to site N + 1 at almost each time step, which leads to ot e tf « 1. Thus, it is impossible for the main chain to maintain 
the low density phase. 

From the analysis above, one can see that there are five possible phases ((LD, LD), (LD, HD), (LD, MC), (HD, HD) 
and (HD, MC)) in this system. Fig. 2(a) shows the possible phase boundaries (a = (3/[m + (m — 1)0\) for m = 1,2 
and 3. With the increase of m, we can predict that the phase boundary moves toward the left in the phase diagram, 
which means that the low-density area decreases while the high-density area increases. The phase diagram for m = 2 
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FIG. 2: (a) Phase boundaries (or transition lines) for m = 1, 2 and 3 in the synchronous TASEPs with an MISO junctions, (b) 
Phase diagram for m = 2 in the synchronous TASEPs with an MISO junction. The solid line represents the (LD, HD) phase 
specified by a = f3/[m + (rn — l)/3] and f3 < 1; the grey oval corresponds to the (LD, MC) phase specified by a = l/(2m — 1) 
and (3 = 1; and the black rectangle for the (HD, MC) phase specified by a > l/(2m — 1) and f3 = 1. 



is also shown in Fig. 2(b). From Fig. 2(b), one can see that: (i) The (LD, HD) phase corresponds to the transition 
line between the (LD, LD) phase and the (HD, HD) phase specified by a — [3/[m + (to — l)/3] and < (i < 1. (ii) The 
(LD, MC) phase is the transition phase between the (LD, LD),(LD, HD), (HD, HD) and (HD, MC) phases specified 
by a — 1/(2to — 1) and (3 = 1 . Also, note that the transition from the (LD, LD) phase to the (LD, HD) phase, the 
density change on the subchains is continuous, while the density change on the main chain is discontinuous. Similarly, 
the transition from the (LD, HD) to the (HD, HD) phases, the density change on the subchains is discontinuous, 
while the density change of the main chain is continuous. Also, for the transition from the (LD, MC) phase to the 
(HD, MC) phase, the density change on the subchains is discontinuous, while the density profile on the main chain is 
unchanged. 

From the above analysis, we can see the non-equilibrium stationary states, stationary-state phases and the phase 
boundary are determined by the boundary conditions of the system as well as by the number of subchains. In other 
words, the dynamics of the system is determined by its environment and its own structure. 
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FIG. 3: Schematic diagram of the domain wall dynamics in the (LD, HD) phase. The domain wall moves in the left and right 
subsystems at rates vl and vr, respectively. 
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III. DOMAIN WALL DYNAMICS 



A phcnomcnological domain wall (i.e., shock) theory to explain phase behavior of a TASEP in a random update 
procedure with open boundaries is introduced in Ref. [36| . A domain wall is a phase boundary connected by two 
possible stationary states. The wall can have a random walk through the system with a drift speed defined as follows 

V= J -±^-. (25) 
P+ - P- 

where J and p are the currents and densities in the two phases; '+'('-') denotes the phase to the right (left) of the 
domain wall. When V > 0, the domain wall moves to the right, while the domain wall travels to the left when V < 0. 
For instance, if a < 3 in a TASEP, the domain wall first appears at the left end and will drift to the right later. The 
wall will pass through the system, which eventually leads to the system being in a low-density stationary state. When 
V = 0, it implies that the domain wall has no net drift between two possible phases. In this case, stationary density 
profiles are linearly increased and the domain wall can exist anywhere with equal probability in the system. 

In this paper, the line specified by a = f3/[m + (m — 1)0\ corresponds to the coexistence of the (LD, LD) and (HD, 
HD) phases. However, the simple approximation theory may not fully reflect the correlations near the junctions as 
indicated in Ref. |23j. The domain wall theory is adopted in order to derive the phase boundary and the density 
profile in the bulk. This theory has also been used in Ref. [231 ] . 

To determine the position of the domain wall in the system, we define x as x — i/(2N), where i is the site index 
and 2N is the length of the system. In the range of < x < 1, the domain wall moves at rate vl in the left subsystem 
(the sub-chains). In the range of 1 < x < 2, the domain wall moves at rate vr in the right subsystem (the main 
chain), see Fig. (3). 



vl and vr can be given by utilizing Eq. ([25 



vl = — J7 j VR = — ^. (26) 

P+ ~9- 9+~ P 



where: 



" L 1 L j §_ R 1 R ' r) -. 

1 + q 1 + a 1 + a 1 + p 1 + p 1+p 



As a result, vl and vr are rewritten as: 



a 3 , , 

W£ = i . v R =-^-r. 28 

1 — a 1 — p 



Similarly to Ref. [23] , (qn) is denoted as a probability to find the domain wall at any position in the left (right) 
subsystem. For a special site, i, in the left (right) subsystem, the probability is obviously equal to qh/N (qn/N). 
Then, at the junction point, we have: 

VlQl = v rQr (2Q s 

In addition, normalized q^ and qn are satisfied with: 

q L + q R = l. (30) 



Combining Eqs. (|29|) and ([30]) . we obtain: 



?L = ; , qn = ; ■ (31) 

VL + VR VL + VR 

According to Eq. Ip8|). Eq. ([31]) implies that: 

^(1-a) _ a(l-8) 

qL ~a + P-2aP' qR ~a + p-2al3- 1 ' 

Accordingly, the probabilities of the domain walls falling in certain zones in the left and right subsystems are also 
given by: 

Prob(xL>w > x) — qLX, < x < 1, (33) 
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and 

Prob(x DW <x) =q L + q R (x - 1), 1 < x < 2. (34) 

Thus, the density at any position in the system becomes: 

p(x)k = P-Prob(xDW > x) + p\Prob(xDW < x), k — L,R (35) 

Finally from Eqs. (|32|) - (f35]) . one can obtain: 

/ s ol 8(1 - a) 2 „ .„„. 

p(x) L = z + T , T7 J n ^ x, < x < 1 (36) 

1 + a (1 + a)(a + 8 - 2<x8) 

and 

n 8 8(l-a)(l~8) a (\-8f 

= TT]3 (1 + 8)(a + 8 - 2a8) + (l + 8)(a + 8- 2a8) {x ~ 1} ' 1 < 1 " 2 ' (3?) 

Densities in the boundary conditions can be calculated as p(a: = 0)l = a/(l + a) and p(x = 2) a , = 1/(1 + 8). These 
results are completely identical with theoretical analysis in Refs. [37l l38l] . At the junction point N, the densities are 
equal to p(x = 1) L = [a 2 (l- 8) + 8(1- a)]/[(l + a)(a + 8-2af3)], p(x = 1) R = 8(1- ad) /[(1 + 8)(a + 8 -2a8)]. Note 
that in the transition line between the (LD, LD) and (HD, HD) phases, we obtain the relationship a — 8/[m+(m— 1)/?]. 



IV. SIMULATION RESULTS AND DISCUSSIONS 



To validate our theoretical analysis, computer simulations are conducted. Here, we only present a synchronous 
TASEP with a Y-type junction, that is m = 2. The numbers of sites of the subchains and the main chain are all 
equal to 1,000. In simulations, stationary density profiles are obtained by averaging 10 8 sampling at each site. The 
first 10 5 iV time steps are discarded to let the transient time out. 

The density profiles for the (LD, LD), (HD, HD) and (HD, MC) phases are shown in Fig. 4. We only illustrate the 
density properties of subchain 1 and the main chain since the density properties of the other subchain is essentially 
the same as subchain 1. It is found that there is a good agreement between Monte Carlo simulations (MCS) and 
mean field (MF) analysis (see Figs. 4(a)-(e)), which verifies our theoretical investigations. 

A phenomenological domain wall (DW) theory developed in Section III is used to calculate the density profiles 
of phase boundaries such as the (LD, HD) and (LD, MC) phases (see Fig. 5). The results obtained from the 
domain wall theory show an agreement with computer simulations. When a and 8 both increase and also maintain 
a = 8/[m + (m — 1)0], the system keeps in the (LD, HD) phase until 8=1; the slope of the density profiles for 
x < 1 decreases until the slope reduces to 0.5, while the slope of the density profiles for 1 < x < 2 also decrease until 
the slope decreases to 0. For instance, the slope decreases from 0.588 to 0.542 (also see Eq. (32)) when a increases 
from 0.1 to 0.2 (see Figs. 5(a) and (b)). Finally, the slopes of density profiles of the subchains become 0.5 and the 
slope of density profile of the main chain becomes when a — 1/3 and 8 = 1 (see Fig. 5(c)). Additionally, Monte 
Carlo simulations, theoretical calculations and domain wall theory all show that, when a = 1/3 and 8=1 (i.e., the 
transition phase between the other four phases), the main chain is in the maximal current phase. 

Density profiles of the systems for m = 2 and m = 3 with the synchronous update scheme are simulated and 
compared (see Fig- E])- According to Eq. (7), the phase boundary between the (LD, LD) and (HD, HD) phases can 
be described as a = 8/(2 + 8) for m = 2 and a = /3/(3 + 2/3) for m = 3. Fig. [S] shows that both systems are in the 
(LD, LD) phase when a = 0.1 and 8 = 0.8. However, when a increases (e.g., a = 0.2), the system for to = 2 is still 
in the (LD, LD) phase, while the system for m = 3 is in the (HD, HD) phase (see Fig. [UJb)). This is due to the phase 
boundary between the (LD, LD) and (HD, HD) phases moving towards the left when to increases (see Fig. 2(a)). 
Density profiles in the (HD, HD) phase for both m = 2 and m = 3 are shown in Fig. HJc). Compared with Fig. [Ha) 
and (c) , it can be seen that the density profiles of the subchains of both systems are the same when both systems in 
the (LD, LD) phase, while the density profiles of the main chains of both systems are the same when both systems 
are in the (HD, HD) phase. Fig. Eld) illustrates that the system is in the (LD, LD) phase for to = 2, while it is in 
the (HD, MC) phase for to = 3. 

We can also see the similarities and differences between the phase diagram of the system with the syn- 
chronous/parallel update scheme (see Fig. 2(b)) and that of the system with the random update scheme (see Fig. 
3 in [23j]). One can see that the structures of the phase diagrams are similar. All have five phases in their phase 
diagrams. Also, increasing the number of subchains (i.e., inputs) only shifts the transition line between the (LD, LD) 
phase and the (HD, HD) phase that does not fall on the boundaries of the phase diagrams. However, the differences 
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FIG. 4: Density profiles obtained from our theoretical calculations and Monte Carlo simulations (MCS) when m = 2: (a) and 
(b) are for the (LD, LD) phase, (c) and (d) are for the (HD, HD) phase, and (e) for the (HD, MC) phase. The parameters are 
set to: (a) a = 0.1 and f3 = 0.8, (b) a = 0.2 and f3 = 0.8, (c) a = 0.8 and /3 = 0.32, (d) a = 0.8 and (3 = 0.44, and (e) a = 0.6 
and = 1.0. 



in the phase diagrams include: (i) the (HD, MC) phase region in the phase diagram of the system with the random 
update scheme reduces to a line in that of the system with the synchronous update scheme; and (ii) the line of the 
(LD, MC) phase in the phase diagram of the system with the random update scheme reduces to a point in the phase 
diagram of the system with the synchronous update scheme. 

Fig. [7] shows the differences in the density profiles of the systems with the synchronous update scheme and the 
system with the random update scheme when m — 2. In Fig. [TJa), these two systems are in the (LD, LD) phase when 
a = 0.1 and j3 = 0.8. When a is increased to 0.2 and (3 is unchanged, the system with the synchronous update scheme 
is still in the (LD, LD) phase, while the phase of the system with the random update scheme becomes the (HD, MC) 
phase (see Fig.JTJb)). Fig.[7^c) shows the system in the (HD, HD) phase in both systems when a = 0.8 and j3 = 0.32. 
With the increase of (3 (e.g., j3 — 0.8), the phase of the system with the random update scheme changes to the (HD, 
MC) phase, while it still keeps in the (HD, HD) phase in the other system. 

Note that the system also exhibits a particle-hole symmetry. Since particles moving forward at junction points 
with the same priority is equivalent to holes moving backward at the same priority. Also, our method can be used to 
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FIG. 5: Density profiles obtained by the domain wall (DW) theory and Monte Carlo simulations (MCS) when m = 2: (a) and 
(b) for the phase coexistence line between the (LD, LD) and (HD, HD) phases, (c) is for the coexistence phases between the 
(LD, LD), (LD, HD),(HD,HD)and (HD, MC) phases. The parameters are: (a) a = 0.1 and (3 = 0.222, (b) a = 0.2 and (3 = 0.5, 
and (c) a = 1/3 and j3 = 1.0. 



analyze synchronous TASEPs with a single-input-multi-output (SIMO) junction. Other inhomogeneous synchronous 
TASEPs can be investigated in the similar way. For instance, it would be interesting to study an MISO junction 
where these parallel subchains are dynamically different. 



V. SUMMARY AND CONCLUSIONS 



Multi-input-single-output (MISO) junctions are relevant to many biological processes as well as vehicular and 
pedestrian traffic flow. Synchronous totally asymmetric exclusion processes (TASEPs) with an MISO junction arc 
investigated in this paper. The theoretical solutions, mean-field approximation, domain wall theory are developed. 
Extensive computer simulations are conducted. Our theoretical analysis suggests that there are five possible stationary 
phases ((LD, LD), (LD, HD), (LD, MC), (HD, HD) and (HD, MC)). The (LD, HD) phase corresponds to the transition 
line (when a = (3/[m + (to — l)/3] and /3 < 1, where to is the number of subchains) between the (LD, LD) phase and 
the (HD,HD) phase. The (LD, MC) phase (when a — l/(2m — 1) and j3 = 1) is the transition phase between the 
(LD, LD), (LD, HD), (HD, HD) and (HD, MC) phases. Also, the non-equilibrium stationary state, stationary-state 
phases and the phase boundary are determined by the boundary conditions of the system as well as the number of 
subchains. The phase boundary moves to the left in the phase diagram when the number of subchains increases. The 
density profiles arc simulated, which shows good agreement with theoretical analysis. 

We also compare the phase diagrams and density profiles between the system with the synchronous update scheme 
and the system with the random update scheme. The main differences in the phase diagrams include: (i) the (HD, 
MC) phase region in the phase diagram of the system with the random update scheme reduces to a line in that of 
the system with the synchronous update scheme; and (ii) the line of the (LD, MC) phase in the phase diagram of the 
system with the random update scheme reduces to a point in the phase diagram of the system with the synchronous 
update scheme. 

The approach used in this paper can be used directly to analyze TASEPs with a single- input-multi-output (MIMO) 
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FIG. 6: Density profiles in Monte Carlo simulations when m = 2 and m = 3. The parameters are: (a) a — 0.1 and /3 = 0.8, 
(b) a = 0.2 and (3 = 0.8, (c) a = 0.8 and f3 = 0.32, and (d) a = 0.25 and (3 = 1.0. 



junction in a parallel updating procedure. 
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